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We study first order phase transitions in the gravitational 
collapse of spherically symmetric skyrmions. Static sphaleron 
solutions are shown to play the role of critical solutions sep- 
arating black-hole spacetimes from no-black-hole spacetimes. 
In particular, we find a new type of first order phase transition 
where subcritical data do not disperse but evolve towards a 
static regular stable solution. We also demonstrate explicitly 
that the near-critical solutions depart from the intermediate 
asymptotic regime along the unstable manifold of the critical 
solution. 

04.25.Dm, 04.40.-b, 04.70.Bw 



The understanding of the dynamics of gravitational 
collapse is a major theoretical challange in classical gen- 
eral relativity. A few years ago a new twist into this 
problem was provided by Choptuik's numerical discov- 
ery of critical behaviour in the gravitational collapse of 
spherically symmetric massless scalar field Q . Analysing 
the evolution of one-parameter families of regular initial 
data, Choptuik showed that the transition between col- 
lapse and dispersion has a universal character which is 
reminiscent of a second order phase transition with the 
continuous change of black- hole mass. Subsequently, sim- 
ilar critical behaviour has been observed in several other 
models of gravitational collapse, and significant progress 
has been made towards the analytical understanding of 
the numerical phenomenology (see for a recent review 
and bibliography). 

Recently, a qualitatively new type of critical col- 
lapse was found in the Einstein- Yang-Mills system ||, 
where, for some initial data, the fundamental Bartnik- 
McKinnon sphaleron (static soliton with one unstable 
mode) plays the role of a critical solution. In this case 
there is a gap in the spectrum of black-hole masses which 
is reminiscent of the first order phase transition. 

Surprisingly, none of the models studied so far had 
a stable static solution (a "star") and consequently the 
possible endpoints of evolution were basically reduced to 
the collapse/dispersion alternative. As far as we know, 
the only exception is the collapse/oscillaton alternative 
recently found in the collapse of massive scalar field jjj. 
If a model has a star-like solution, then there will be an 
open set of initial data which evolve towards this solu- 
tion (basin of attraction) and it is natural to expect the 
occurence of critical behaviour at the boundary of this 



set. 

The Einstein-Skyrme model which possesses the whole 
zoo of static regular solutions, both stable and unstable, 
seems to be a suitable testing ground to address this is- 
sue. In this letter we focus our attention on a new type 
of first order phase transition we found in this model for 
baryon-number-one data. We also briefly discuss a novel 
feature of mass gap in the baryon-number-zero sector. A 
detailed description of both first and second order phase 
transitions will be published elsewhere 

We consider the Einstein-Skyrme system, so the mat- 
ter in our model is an SU (2)-ve\\iec\ scalar function U(x) 
(called a chiral field) with dynamics given by the La- 
grangian 
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L = J — Tr(V a V a U~ 1 ) + —!-^Tr[(S7 a U)U~ 1 ,(V b U)U~ 
4 326^ 

where V a is the covariant derivative with respect to the 
spacetime metric. The two coupling constants / and e 
are hereafter set to one, which amounts to using f/e and 
1/e/ as the units of mass and length, respectively. 

We specialize to spherical symmetry. For the metric we 
use the polar time slicing and the areal radial coordinate 



ds 2 



. e -2S(r,t) N ^ ^ dt 2 + N -i( r , t)dr 2 + rW. (2) 



For the chiral field we assume the hedgehog ansatz 
U = exp(ia ■ fF(r,t)), where a is the vector of Pauli 
matrices. Using overdots and primes to denote d/dt 
and d/dr respectively, we introduce an auxilary variable 



P = ue d N 1 F, where u 



2sm z F. Then, the 



Einstein-Skyrme equations reduce to 

F = e~ s N—, 
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Here a = AirGf 2 is the dimensionless coupling constant. 
In order to solve the initial value problem numerically wc 
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have implemented a free evolution scheme using a sec- 
ond order finite difference method on uniform grid. The 
conservation of the hamiltonian constraint (7) was used 
only to check the accuracy of the code. To ensure reg- 
ularity at the center we impose the boundary condition 
F(r,t) = 0{r) for r — » 0. To minimize reflections from 
the outer boundary of the grid we impose a sort of out- 
going wave condition there. Asymptotic flatness requires 
that F(r, 0) = Bir + 0(l/r 2 ) at infinity, where the integer 
B, called the baryon number, is equal to the topological 
degree of the chiral field. As long as no horizon forms, the 
baryon number of initial data is automatically preserved 
during the evolution. 

Below we describe our results for regular initial data 
with baryon number one and zero. 

B = 1 : Let us recall || that in this case, for a < a% ~ 
0.040378, the Eqs.(3-7) have two regular static solutions 
(solitons) which we denote below by X s and X u with X 
standing for (F, P, N, 5). The solution X s is linearly sta- 
ble - it is the skyrmion distorted by gravity. The second 
solution X u has one linearly unstable mode. These two 
solutions coalesce at the bifurcation point a — a±, and 
disappear for a > a±. For given a the mass m u of the 
unstable soliton is larger than the mass m s of the stable 
soliton (for example, for a = 0.02 we have m s ~ 68.05 
and m u ~ 80.70). As a — ► 0, m s tends (from below) to 
the mass of the flat-space skyrmion (~ 72.92), whereas 
m u — » oo, which is the manifestation of a nonperturba- 
tive character of the solution X u . 

Since for B = 1 the dispersion of the chiral field to 
infinity is forbidden, generic initial data end up either as 
black holes or as the stable soliton X s . The respective 
basins of attraction of these two final states depend on 
the coupling constant a. This is illustrated in Fig. 1 for 
a typical time-symmetric kink-type initial data. 
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FIG. 1. Basins of attraction of the one-parameter family 
of B — 1 initial data F p (r,0) — 7rtanh(r/p) as a function of 
the coupling constant a. Critical initial data are represented 
by the solid line curve. Initial data lying in the region en- 
closed by the critical curve and the p-axis evolve towards the 
skyrmion X s , while initial data lying outside that region form 
a black hole. As a — > (not shown in the figure) the basin of 
attraction of X s tends to the whole positive p-axis. 



Notice that for a < a± there are two critical config- 
urations along the one-parameter family F p (the points 
A and B on the dashed line segment on Fig. 1). Since 
the parameter p (the width of a kink) is inversely pro- 
portional to the "condensation" of initial mass, the 
transition at the point A could be interpreted as a 
"weak/strong" field transition. However at the point B 
it is the "weaker" configuration that collapses, simply be- 
cause it is unable to get rid of the excess mass by the time 
it shrinks to size of the stable soliton. This behaviour 
indicates that it is difficult to formulate a criterion for 
collapse in terms of initial data alone. 

The critical initial data lying on the boundary of the 
basin of attraction of X s asymptote to the unstable so- 
lution X u . In other words, the codimension one stable 
manifold of the solution X u divides (locally) the phase 
space into collapsing and non-collapsing data. Near- 
critical data approach X u , stay in its vicinity for some 
time, and eventually collapse to a black hole or decay into 
X s . This is shown in Fig. 2 for a marginally subcritical 
data (lying slightly below the point B on Fig. 1). 
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FIG. 2. A sequence of profiles of F(r, t) (solid lines) versus 
fn(r) from a marginally subcritical evolution of the B — 1 
initial data. The stable and unstable static solutions are rep- 
resented by the dotted and dashed lines, respectively. Here t 
is the proper time at spatial infinity (we use the normalization 
<S(oo,t) = 0). 

During the final stage of the evolution of a marginally 
supercritical solution (i.e. when the solution runs away 
from the vicinity of the critical solution and forms a 
black hole) almost no energy is radiated away to infinity. 
Therefore the mass gap in this first order phase transition 
equals m u (up to 0.1%). 

For subcritical data the process of settling down to 
the stable soliton seems to have a universal late-time be- 
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haviour which is dominated by a fundamental quasinor- 
mal mode. The parameters of this mode depend strongly 
on a. Having convinced ourselves that these parameters 
are indeed universal, that is do not depend on the ex- 
citation, we have determined them from the evolution 
of finite perturbations of X s (see Fig. 3). The details 
of such perturbation are washed out very rapidly and a 
characteristic exponentially damped " ringing" dominates 
the late time evolution. 
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FIG. 3. Quasinormal ringing about the stable soliton for 
a = 0.03. We plot dN = \N(r ,t)-N s (r )\ from the evolution 
of initial data F(r, 0) = 7rtanh(r) (dots). Here ro = 0.8. The 
least squares fit for the times t > 10 (represented by the solid 
line) gives the frequency uj ~ 0.494 and the damping rate 
r ~ 11.69. 

The above results are in perfect agreement with lin- 
ear |?],|| and nonlinear || stability analysis of gravitat- 
ing skyrmions. In the intermediate asymptotics the evo- 
lution of near-critical data is well approximated by the 
linearization about X u 



X p (r,t) « X u (r) + C(p-p*)e xt 6X u {r) + 



decaying 
modes ' 



(8) 



where 8X u (r) is the single unstable eigenmode associ- 
ated with a positive eigenvalue A. Depending on whether 
p > p* or p < p*, the solution X p eventually collapses 
to a black hole or decays to X s . The "lifetime" T of 
a near-critical solution staying in the vicinity of X u is 
determined by the time in which the amplitude of the 
unstable mode grows to a finite size: \p — p*\e XT ~ 0(1), 
which gives T ~ — A -1 In \p — p*\. Thus, the larger A, 
the better fine-tuning is required to see the solution X u 
clearly pronounced as the intermediate attractor. 

We have verified the formula (8) by comparing the 
unstable eigenmode 6X U with the snapshots of a near- 
critical solution departing from the vicinity of X u . This 
is shown in Fig. 4 for the same initial data as in Fig. 2. 
For t — 24 the profile of F(r, 24) is practically indistin- 
guishable from the (suitably normalized) unstable eigen- 
mode SF u (r). A slight deviation of F(r, 23) from SF u (r) 
for large r is due to the fact that by t = 23 the decaying 
modes in (8) have not yet died away completely. On the 



other hand, for t = 25 the deviation for small r signals 
the onset of nonlinear regime. 
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FIG. 4. The snapshots of F(r, i) (solid lines) departing 
from the vicinity of the critical solution X u . On each snap- 
shot we superimpose (dashed lines) the unstable eigenmode 
5F u (r) computed via linear perturbation analysis. For t — 24 
the amplitude of 8F u (r) is normalized to the maximum of 
F(r, 24). The corresponding amplitudes for t = 23 and t — 25 
are rescaled by the factors e _A and e x , respectively. 

Using (8), we have also computed the eigen- 
value A directly from the nonlinear evolution by 
monitoring F p (r, 24) at several discrete radii r^, 
i = l,...,n, and evaluating the averaged quantity 
J2 In (Ppin, 24 + At)/F p (n, 24)) for some small At. 
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The result agrees up to three decimal places with the 
eigenvalue obtained via linear stability analysis. 

B = 0: In this topological sector, the Eqs.(3-7) have a 
pair of static regular solutions for a < ao ~ 0.00147 ||. 
One of them has two unstable modes, hence it plays no 
role in the evolution of generic initial data. The sec- 
ond solution, call it Y u , has one unstable mode, so it is 
a candidate for the critical solution. Indeed, analysing 
the ingoing "generalized Gaussian" profile F(r, 0) = 
^ r 3g-(r-r (l ) /A ^ we j-^yg f ounc [ y u as the intermedi- 
ate attractor at the border between collapse and disper- 
sion. Since in this case the critical behaviour is anal- 
ogous to the first order phase transition in the EYM 
model here we discuss only the issue of mass gap in 
this transition. Analysing near-supercritical solutions we 
found that, in contrast to the previously studied cases, 
the mass gap is not equal to the mass of the sphaleron. 
In order to determine the mass gap more precisely we 
adopted the following strategy. Rather than improving 
the accuracy of fine-tuning (which is computationally 
time-consuming because the eigenvalue of the unstable 
eigenmode is rather large), we began the evolution with 
specially prepared initial data of the form of the static 
solution Y u plus the unstable eigenmode with a small 
positive amplitude: Y(r, 0) = Y u (r) + e8Y u (r). In other 
words, we pushed the static solution along its unstable 
manifold and let it collapse. It turned out that the mass 
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gap is considerably smaller than the mass of the critical 
solution (which of course sets the upper bound for a mass 
gap). This is illustrated in Fig. 5, where one can clearly 
see how a substantial amount of the initial mass is being 
radiated away to infinity. The resulting black hole has 
the mass ~ 253.1, which is 75.6% of the mass of Y u . 
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FIG. 5. The mass aspect fi(r, t) (defined by TV = 1 — 2/i/r) 
from the evolution of the static solution Y u perturbed along 
the unstable direction (for a = 0.00145). The mass inside 
radius r at time t measured in units of f/e is given by 
m(r,t) — —fj,(r,i). The plateau of fJ.(r,t) determines the 
mass gap. 

Concluding, we believe that the results presented here 
have clarified certain aspects of first order phase transi- 
tions in gravitational collapse and, in combination with 
the results of |3| and Q, they help to understand which 
features of this phenomenon are generic. 
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